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ABSTRACT 

Oceanographic and meteorological data from Ocean Weather 
Station (OWS) NOVEMBER were analyzed by statistical and 
spectral methods in order to describe the nature and perio- 
dicities of the air-ocean energy exchange. Salinity data 
from 1968 through 1970 were compared with daily observed 
parameters commonly associated with changes in the salinity. 
Regression analysis showed surface salinity to be highly 
correlated with surface pressure which is interpreted as the 
surface salinity responding to a baroclinic transport. At 
this location none of the classical parameters usually 
associated with fluctuations in the surface salinity were 
highly correlated. 

Twenty-four years (1947-1970) of daily averaged, 
observed climatological values, as well as heat exchange 
terms computed from them, were compared with the sea-surface 
temperature time series and spectrally analyzed. The 
spectral analysis yielded spectral density, phase and 
coherence functions indicating the response of surface 
temperature to the heat exchange terms. From these it was 
evident that at the yearly cycle the net heat exchange, whi ch 
is shown to be equal to the horizontal heat advection, is 
of the same order of importance as solar radiation in its 
effect on the fluctuation of sea-surface temperature. 
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I. INTRODUCTION 



A. OBJECTIVE 

Sea-surface temperature is one of the most important 
parameters in any study of air-sea interaction. Anomalies 
of this parameter, measured over periods of a week to 
several months, have been linked to gross atmospheric 
effects. However, few studies have investigated the time 
response of the sea-surface temperature over long periods 
to those factors which govern its fluctuations. 

Until the more recent decades such a study would have 
been limited by both instrumentation and general paucity 
of data. This resulted in reliance upon monthly or annual 
means to describe the state of the ocean. The advent of the 
ocean weather station program in the 1940's, and recently, 
the use of satellites has greatly increased the amount and 
frequency of data available to the researcher. Also, the 
advent of high speed computers has made the use of time 
series analysis routines relatively easy and inexpensive to 
use . 

The past limitations made statistical studies of changes 
of surface temperature on time scales of a day or less, 
difficult. However, daily fluctuations can be quite inform- 
ative. If advection is ignored, the sea-surface temperature 
at OWS NOVEMBER was found to fluctuate primarily in response 
to the day to day changes in air temperature, the amount of 
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insolation, wind-wave mixing, and the amount of cloud cover 
as well as the elevation of the sun. The magnitude of 
these diurnal fluctuations can be as large as the monthly 
f 1 uctuations . 

Sea-surface temperature fluctuations are important in 
many oceanographic parameters, such as the stability of the 
water column. This, in turn, is of great consequence in 
sound propagation. Anomalies of sea-surface temperature, 
likewise, have been linked to climatic variations over the 
ocean . 

The salinity is important in that it affects the 
density of sea water. Thus, it affects both vertical 
stability and sound propagation. Together with the water 
temperature, one can infer transport characteristics of a 
region based on water mass analysis. 

This study examines changes occurring at the surface 
of the ocean over time scales varying from one day to 
several years. The two parameters of primary interest are 
the sea-surface temperature and surface salinity. The area 
of study is Ocean Weather Station (OWS) NOVEMBER located at 
30° North and 140° West. The data used are synoptic data 
collected by the United States Coast Guard at this position. 
The data are of two types: hydrographic Nansen cast data and 
shipboard meteorological observations. 

The thesis is organized to first present correlations 
of the primary parameters versus various meteorological 
parameters. These correlations provide an indication of 



those parameters which can be used as predictors for the 
sea-surface temperature and salinity. 

Secondly, a spectral analysis of the data is made. 

From this, the synoptic atmospheric mechanisms which cause 
the fluctuations of the sea-surface temperature are identi- 
fied. Periods of high spectral energy, along with their 
corresponding phase and coherence relationships are examined 
to gain insight into the time scales associated with each 
mechanism. 

B. REVIEW OF THE LITERATURE 

1 . The Subtropic Region 

Ocean Weather Station NOVEMBER, located at 30° 

North and 140° West, lies within what Tully (1964) defines 
as the Pacific Subtropical Region. Figure 1 shows its 
location in comparison to other major features of the North 
Pacific Ocean. The region is characterized by three dominant 
conditions: an excess of evaporation over precipitation 
throughout the year, the occurrence of convective mixing, 
and seasonal heating and cooling patterns which result in a 
pronounced seasonal thermocline. 

The excess of evaporation over precipitation (E-P) , 
has a distinct effect on the other two traits. Tully (1964) 
estimated E-P to be approximately 80 cm” ^ yr in the vicinity 
of OWS NOVEMBER. This leads to a removal of fresh water 
from the surface and acts as a driving force for the con- 
vective mixing. This mixing then contributes to the erosion 
of the seasonal thermocline. 
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2 . 



Surface Circulation 



The subtropical region of the Pacific is dominated 
by the Pacific Gyre. Figures 2 and 3 show graphically how 
the gyre varies seasonally both in location and intensity. 

In the winter the gyre's axis is displaced northward; in 
summer it is shifted approximately 100 n mi southward and 
intensifies by approximately 0.2 kt. Throughout the year 
the current speed ranges from 0.1 to 0.3 kt. Hansen (1973) 
hypothesized that this intensifi cation occasionally results 
in the advection of large volumes of cold subarctic water 
into the area surrounding 0WS NOVEMBER, specifically during 
the period of 1968 through 1970. 

3 . The Subarctic Transition Zone 

The North Pacific Gyre, linked with the Subarctic 
Gyre, is also involved in the formation of a transition zone 
which bounds the subtropical region to the north. Roden 
(1970) first described the subarctic transition zone as an 
area characterized by large temperature and salinity gradients 
and low hydrostatic stabilities. This oceanic front defines 
the boundary between the subtropical and subarctic water 
masses, varying in location between 30° and 40° North. 
Subsequent articles by Roden (1971, 1972) further described 
this zone in the North Pacific. 

The transition zone is largely maintained by the 
surface wind stress. Its northern boundary is characterized 
by the subarctic water mass, resulting in temperature 
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Figure 2. Winter surface circulation patterns for the North 
Pacific Ocean (Department of Commerce, 1961). 
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Figure 3. Summer surface circulation patterns for the North 
Pacific Ocean (Department of Commerce, 1961). 



inversions, gravitational instabilities, a complicated sound 
velocity structure, and a disappearance of the subarctic 
halocline. It appears to lie in the area of 40N to 45N 
latitude throughout the year. 

The southern boundary lies in the area of 30N to 
32N, depending upon the season, and is displaced to the 
north more in the winter months. It is characterized by 
sharp temperature gradients in the thermocline and vertical 
instabilities in the mixed, upper layers. A salinity minimum 
is a common occurrence at this boundary. Also, around 30N 
there is a transition toward the higher salinities and 
higher temperatures of the subtropical water mass with 
salinities greater than 35 0 / oo and a strong halocline located 
at a depth between 50 m and 100 m. 

Thus, the transition zone is defined by a deep 
salinity minimum, at approximately 500 m and a concurrent 
sound velocity minimum of less than 1480 m sec"^. There is 
apparently both a seasonal and yearly variation in the posi- 
tion of the boundary (Seckel, 1968). Both the salinity 
minimum and sound velocity minimum rise as one goes further 
to the north. Roden (1970) also examined both barotropic 
and baroclinic motion in the transition zone and found that 
they rarely exceeded 5 cm sec - ^. 

As the transition zone approaches North America, 
it turns to the southeast and becomes associated with the 
heat and salt flux divergence zones noted by Roden (1971) 
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in the area of OWS NOVEMBER. These are related to flux 



divergencies which are linked to the presence of horizontal 
advection of salinity as has been noted by Hansen (1973). 

4 . Average Annual Heat Balance 

Based upon an examination of the annual heat 
balance of the North Pacific Ocean and the prevalent circu- 
lation patterns, Wyrtki (1965) determined that there is a 
small amount of heat gain in the vicinity of OWS NOVEMBER. 

On the average, he estimated this gain to be approximately 
25 ly day (where a langley, ly, = 1 cal cm ). 

Seckel (1970) computed the heat balance for the 
area of OWS NOVEMBER during the Trade Wind Zone Pilot Study. 
His studies showed that the average net radiation computed 
on a daily basis was approximately 3 ly day"^. 

In a more recent study, Dorman (1974) studied 20 
years of the heat balance at OWS NOVEMBER and separated it 
into two decades, 1951 through 1960 and 1961 through 1970. 

He found the two decades had distinctly separable heat 
balances, with the first decade having a deficit of 
90 ly day"^ and the second a 40 ly day"^ deficit. These 
results were based on three-hourly meteorological observa- 
tions. Dorman ascribed this difference to the change in 
climatic regimes described by Namias (1970, 1971), which 
resulted in displacement of cooler waters in the North 
Central Pacific by warmer water, beginning approximately in 
1961 . 
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5. 



Surface and Subsurface Temperature Anomalies 



The sea-surface temperature and its anomalies have 
been widely studied in their relationships to both oceano- 
graphic and meteorological forecasting. Emphasis has also 
been placed on subsurface temperature anomalies. Beland 
(1971) illustrated their considerable importance to the flow 
of energy between the sea and the atmosphere. Bernstein 
(1974) used subsurface temperature measurements in the trade 
wind region northeast of Hawaii to identify mesoscale eddies 
drifting westward. 

Namias (1971) in his study of an anomalous warming 
of the northeastern Pacific from May to June 1968, found the 
anomaly extended to a depth of 122 m. The warming, as a 
whole, was attributed to increased insolation and horizontal 
convergence of surface waters with its associated down- 
welling. The warming was identified with the development 
and maintenance of a Pacific anticyclone in June of 1968. 

The depth of the anomaly was attributed to bulk heating and 
downward advection of the warmer surface waters. 

Clark (1972) concluded that both advection and 
non-advecti ve processes are vital in determining the nature 
of surface temperature anomalies. He hypothesized that 
advection is the more important process in the winter and 
spring seasons, while the non-advecti ve processes dominate 
in the summer and fall. 
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Beland (1971) in his study of the North Pacific, 
observed that about 50% of the near-surface temperature 
anomalies extended to at least 100 m depth. Most of these 
were negative anomalies in the months of March and April, 
apparently due to the distribution of heat throughout the 
upper layer by convective mixing. Positive anomalies were 
more common in December through February and were confined 
to the upper 30 m. Whereas the negative anomalies were of 
a magnitude of 1.0°C or more, the positive anomalies were 
generally less than 1.0°C. 

Finally, Favorite and McClain (1973) reported that 
they could observe an orderly trans -Paci f i c movement of 
areas of anomalously warm and cold surface waters. These 
features were found to be quasi -permanent and apparently 
move around the North Pacific gyre, approximately 180 degrees 
out of phase and with periods of 5 to 6 years. These 
anomalies have been shown to have a great influence on the 
weather conditions off the west coasts of Canada and the 
United States, as shown by Namias and Born (1970). 

From the above it can be noted that temperature 
anomalies, both surface and subsurface, and whether temporary 
or periodic, inhabit the region around 0WS NOVEMBER, and in 
all probability, affect the location of the subarctic 
transition zone and the heat balance through the process of 
advecti on . 
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6 . 



Recent Studies of Pacific Air-Sea Interactions 



In recent years there have been three studies 
of the air-sea relationships at or around OWS NOVEMBER by 
different researchers using three different methods. Taken 
as a whole, they provide a logical progression, resulting 
in the methods utilized in this study. The first of these 
studies was carried out by Seckel and Yong in 1970, using 
harmonic analysis to study temperature and salinity fluctua- 
tions at Koko Head, Hawaii, and Christmas Island. Hansen 
(1973) used regression methods to study surface and 
subsurface temperatures and salinities at OWS NOVEMBER. 

Dorman (1974) used a spectral approach on time series of 
meteorological values at OWS NOVEMBER, dividing the time 
series into tw,o ten-year decades. Since each of these 
studies utilized a different approach, each will be described 
in some detail here. 

Seckel and Yong (1970) utilized temperatures and 
salinities measured at two points in the North Central 
Pacific; Koko Head, Hawaii, beginning in 1956 and running 
through 1969, and Christmas Island, beginning in 1954 and 
continuing through 1969. Their hypothesis was that 
fluctuations in the temperature and salinity time series 
reflected air-sea interaction processes (heat exchange, E-P) 
and oceanographic processes (advection, diffusion, etc.). 
Their approach used a Fourier series to provide a least 
squares fit to the observed values, which was done for each 
year of the available data. This approach expressed the 
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temperature and salinity as a function of time. Phase 
angles and coef f i ci ents -for each year of salinity and temper 
ature data from both locations were computed and tabulated. 

The results at Koko Head showed seasonal varia- 
bility in the salinity, but not the temperature series. At 
both locations the annual cycle (n=l) was the dominant cycle 
Also of note were the 60- and 52-day cycles, ascribed by the 
authors to the presence of large geostrophic eddies. The 
Christmas Island results were quite similar except that they 
lacked evidence of the eddy-induced seasonal variations. 

The next major method was a regression analysis 
used by Hansen (1973) in his study of temperature and 
salinity fluctuations at the surface, 50 m, and 200 m at 
OWS NOVEMBER. His data were for 1968 through 1570 and 
consisted of daily hydrographic records. Deterministic 
approaches showed that the temperature system was seasonal 
so the data were divided into four seasonal records, 

February to April, etc., in three month segments correspond- 
ing to Spring, Summer, etc. Regressions were then computed 

* 

with the surface salinity, and 50 m and 200 m temperatures 
being regressed against the surface temperatures. 

Most of the results showed little correlation, 
except occasionally during the spring and fall transition 
periods. Only a very poor correlation was found between the 
surface salinity and the surface temperature. From this, it 
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was concluded that these seasonal correlations were due to 
the advection of water by the shifting of the subarctic 
transition zone. 

The final approach was that utilized by Dorman 
(1974) who computed spectra for surface wind speed, air 
pressure, sea temperature, air temperature and dew-point 
temperature for the 10-year period 1961 through 1970. 

Spectra were also computed for each of the four seasons. 

From this it was found that the energy peaked significantly 
at the annual, diurnal and semi-diurnal periods. Variances 
were largest in the winter season and smallest in the 
summer season, except for the sea-surface temperature which 
had its largest variance in the spring and the smallest in 
the winter. The maximum spectral density, for all variables, 
was at the annual period. From the annual period the 
spectral densities followed a decrease to the shorter periods, 
where significant spikes at the diurnal and semi-diurnal 
periods were noted in many of the spectra. 
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II. TREATMENT OF THE DATA 



A. SOURCES 

The climatological data used in this study were taken 
from the NWRC, Asheville data file of ocean weather station 
observations. These observations were begun at OWS NOVEMBER 
in 1947 and have continued to the present. The observations 
were taken every three hours, except for some short periods 
when all or only some of the observations were not taken for 
some parameters. The data extracted included the sea-surface 
temperature, air temperature, wind speed, surface pressure, 
cloud amount and the present weather codes for 1947 through 
1970. 

The salinity data used in this study came from the NODC 
hydrocast data. Since July 1966, the U.S. Coast Guard has 
conducted a continuing program of collecting Nansen cast 
data at OWS NOVEMBER. The surface salinity data for 1968 
through 1970 was utilized in one segment of this study. 

B. TREATMENT OF THE DATA TIME SERIES 

Daily surface salinity values from 1968 though 1970 
were extracted from the NODC data set and any gaps in the 
series were linearly interpolated. The daily salinity values 
were then included with surface pressure, evaporation, solar 
radiation, heat of evaporation, and wind speed, all taken 
directly from or computed from the Asheville data set. To 
this was added the calculated values of precipitation and 
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the evaporation minus precipitation terms to form a three- 
year time series which was then used in the regression 
analysis segment of this study. 

The Asheville meteorological data were used primarily 
to create the 24-year (1947-1970) time series used in the 
spectral analysis segment of this study. The values of SST , 
air temperature, cloud amount, air-sea temperature difference, 
the saturated vapor pressure of the air, the vapor pressure 
over the water, wind speed, surface pressure and wind speed 
squared were extracted or computed for every three hours 
and then converted into daily average values and these were 
the values utilized in all the spectral analyses. For the 
cloud amount parameter, only daylight observations were used. 
From these data sets the heat budget terms were computed on 
a daily basis. For the occasional gaps in the data, linear 
interpolations were made. 

C. THE HEAT BUDGET EQUATIONS 

Empirically derived formulas were used to create the 
heat flux data for this study. The heat budget equation 
utilized excludes the effects of heat gained or lost by 
ocean currents : 



HEAT CHANGE = HEAT INCOMING - HEAT LOSSES 
QN = QS - QE - QB - QC 



( 1 ) 
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where QN is the measured heat change, i.e., the heat gained 
or lost by the water due to the exchange processes. On the 
right-hand side of the equation, QS is the incident solar 
radiation, QB is the effective back radiation, QE is the 
heat of evaporation and QC is the conduction of sensible 
heat. Each of these terms will be described in detail below 
as to their formulation. 

The equation used for QS is an improved version of that 
described by Seckel (1970). The basic form is: 

QS = Q q K (2) 

whe re 

Q q = A q + A-| cosij) + B-j sinif) + k^ cos2<j) + B 2 sin2(j> (3) 

From the Smithsonian Tables (List, 1951), utilizing a 
transmission coefficient of 0.7, the above terms can be 
computed as follows for latitude 30N: 



> 

0 

II 


-32.65 


+ 


674.76 


cos a 


= 


551 . 7 


(4) 


A 1 - 


19.88 


+ 


397.26 


cos 


(2+90) = 


-178.7 


(5) 


B 1 - 


- 6.75 


+ 


224.38 


s i n 


z 


105.4 


(6) 


A 2 


-1.32 


+ 


16.10 


s i n 


2(2-45) = 


- 9.4 


(7) 


B? 


-1.04 


+ 


29.76 


cos 


2(i-5) = 


18.1 


(8) 



where : 

Jt is the latitude 
t is the time of the year in days 

C is the cloudiness in tenths, ranging from 0 to 1 , and 
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(9) 



♦ " 365 ( t ~ 21 ^ 



The cloud factor correction utilized was one for 
cumulus type clouds (Seckel, personal communication): 



An albedo factor of 0.95 was utilized. 

The new equation results in values approximately one- 
fourth to one-third greater than previously estimated. 

The effective back radiation, QB, is a Berlyand 
expression utilized by Seckel (1970): 



where T w is the sea-surface temperature in degrees C. It 
utilizes a cloud coefficient of 0.6 and neglects latitudinal 
variation. 

QE, the heat of evaporation depends upon the difference 
between the saturation vapor pressure over the sea at the 
sea-surface temperature and the vapor pressure of the air 
above the sea surface: 



K = 1.0 - 0.012C - 0.363C 



2 



( 10 ) 



QB = 1.14 (10)' 7 (273.16 + Tj 4 



( 11 ) 



(0.39 - 0.05 V^) (1 - 0.6C 2 ) 



+ 4.58 (10)" 7 (273.16 + Tj 3 




QE = 3767 C n (0.98 e -e ) W 
D w a 



( 12 ) 
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where : 



e is the vapor pressure of the air (mb) 

Cl 

e is the saturation vapor pressure over water 
w at sea-surface temperature (mb) 

Cp is the drag coefficient 

W is the wind speed (m sec"^) 

Whereas Seckel utilized a varying drag coefficient, this 

_ 3 

study utilized a constant drag coefficient of 1.3 x 10 
As in Seckel 's study, the vapor pressure difference was 
multiplied by 0.98 to account for the salt effect. 

The final heat balance term is the conduction of 
sensible heat, QC: 

QC = 4 (T w -T a ) W (13) 

where : 

T w is the sea-surface temperature (C) 

T is the air temperature (C) 

a 

W is the wind speed (m s e c ~ ^ ) 

It can be seen readily that QC is dependent upon two factors, 
the air-sea temperature difference and the wind speed. It 
is usually the smallest of all the heat exchange terms. As 
with Seckel 's study, a conduction factor of 4.0 is used. 

The above equations were those used to create the five 
heat flux time series used in this study. QS was computed 
using only the daylight cloudiness observations, night 
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values being ignored. A daily value for each day was 
calculated for the span of 1947 through 1970. They were 
compared with Seckel's (1970) results to guarantee their 
accuracy, where possible. The utilization of the newer QS 
formula resulted in higher values for this term than he 
reported and correspond!' ngly more positive values of QN . 

D. COMPUTATION OF PRECIPITATION 

In relation to the study of the surface salinity at 
OWS NOVEMBER, it was vital to have some measurement of E-P, 
the e vaporati on-preci pi tati on term. The evaporation was 
readily available, having been calculated from the data 
already available. However, no direct measurement of pre- 
cipitation is made routinely at ocean weather stations. 

To gain a measure of the precipitation, Tucker's (1961) 
method was used. Elliot and Reed (1973) applied this method 
along the Pacific Northwest Coast with good results. This 
approach is based on the current weather code which gives a 
numerical value to the type of weather being observed at 
three hour intervals; this yields an estimation of the amount 
of rainfall. The amount of rainfall was based on the amount 
of precipitation observed for a similar weather code at a 
land station where actual rainfall measurements were avail- 
able. Based on a frequency analysis for all the codes and 
the measured rainfall. Tucker was able to fit a numerical 
system of equations to the prediction of precipitation, 
using only the weather code as an input. 
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In Tucker's system all codes less than 50 are rejected 
since they contribute no appreciable amount of precipitation. 
The codes used had values ranging from 50 to 99. Each of 
these codes is associated with a term made up of some combi- 
nation of three rainfall parameters. These are x, associated 
with light, continuous rainfall; y, associated with moderate, 
continuous rainfall; and z, which is a measure of heavy, 
continuous rain. These three parameters satisfy the follow- 
ing criteria: 

y = nx mm where 1 < n < 5 
z = my mm where l<m<3 
z = 10 mm (in 3 hr) 

For this study x was set equal to 1.85, y = 5.66, and 
z = 8.13, the same values as derived by Tucker from his 
frequency distributions of the codes recorded and the actual 
rainfall at land stations. Thus, the amount of rainfall, 
in millimeters, associated with each current weather code 
can be determined from one of the above, or a combination or 
ratio of these terms. For example, weather code 64, which 
is "rain, not freezing, intermittent, but heavy at the time 
of observation", is described by the term z/2. Likewise, 
code 84, "showers of rain and snow mixed, moderate to heavy", 
is represented by the formula, . 
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For the purposes of this study, the precipitation was 
calculated every three hours and then summed to obtain the 
total precipitation for that particular day. This was then 
subracted from the evaporation term for that day to yield 
the E-P time series for 1968 through 1970. 

E. COMPUTATION OF THE ENERGY SPECTRUM 

The covariance function, R (x), of a time series of 

A A 

random data is a measure of the general dependence of the 
values of the data at one time, t, with the values at 
another time, t+r. The covariance equation in integral form 
i s : 

T 

R xx (t) = lim J x (t) x ( t+x ) dt (14) 

T-*» o 

This yields the average of the product of the value of a 
data point, x(t), with that of x(t+x), a point at some 
lagged time increment, x. As T approaches infinity, the 
resulting average product approaches the exact auto- 
covariance function. 

Since this data forms a finite data set, a lag window 
is applied to the covariance function. The Parzen window 
used in this analysis is defined as: 



P(t) - 1-6( T /T m r + 6(x/T m ) J 


x = 0,1,2 --- T / 2 

m 


(15) 


- (Mt/T „» 3 


x = T /2 + 1 , ---T 
m m 


(16) 


= 0 


x = T m 
m 


(17) 
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where x is the lag and T m is the total lag time. This 
window has the desirable effect of minimizing any numerical 
instability due to side lobe resonance (Bendat and Piersol, 
1971 ) . 

The energy-densi ty spectrum is obtained by Fourier 
transforming the covariance function. The form of the 
energy-density spectrum, corrected by the Parzer. window is: 

oo 

S xx (f) = J P(t) R xx ( t) e " i2nfx <^ (18) 

— 00 

with f being a specific frequency. The energy-density 
spectrum is the statistical estimate of the energy-density 
distribution for the frequency bands forming the total 
spectrum. 

F. COMPUTATION OF THE CROSS-SPECTRAL DENSITY 

When comparing two-time series of the same length and 
time interval, the cross - covari ance function similar to the 
covariance function, is used. For a given lag time interval, 
x , it has the form : 

T 

R xy (f) = lim | | x (t)y (t+x) dt, (19) 

T-5-00 o 

Because the data sets are records of finite length, the 
Parzen lag window is applied. Thus, similar to the covariance 
function and the spectral density, the cross-spectral density 
is the Fourier transform of the Parzen window lagged, cross- 
covariance function. 
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The cross-spectral density is a non-zero function, 
having real and imaginary parts of the form: 



V f) 


00 

= / PM S xy (x) e - f2l,f T dt (20) 

- 00 



which can be expressed in terms of its real and imaginary 



parts : 




V f) 


- C xy ( f ) - IQ <f> (2D 



where C (f), the real co-spectrum, is given by: 
xy 



V f) 


00 

= 2 J [S X y (t ) + S X y(-x)] cos ( 2 7T f T ) dx (22) 

0 



and the imaginary quadrature spectrum, Q (f): 



V f) 


OO 

= 2 f [ S xy ( x ) - S xy ( - x ) ] sin ( 2 tt f x ) dx (23) 

0 



G. COMPUTATION OF THE COHERENCE FUNCTION 

The coherence function, for two records x and y, is a 
measure of the statistical dependence of two series. The 
coherence is defined as: 



2 

Y ,( f ) 
xy v ' 


s ff s (TT 

xx v J yy v ; 


where : 





0 < y 2 (f) < 1 . 

xy 
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are totally coherent and statistically dependent. If 
2 

y (f) = 0 for all frequencies, then the x and y records 

xy 

are totally independent. For intermediate range coherences 
three conclusions may be inferred concerning the two records 
(Bendat and Piersol, 1966); (1) the systems are non-linear, 
(2) noise is present in the measurements, or (3) the output 
response is due to more than one input function. It should 
be noted that this is the coherence squared, but will here- 
after be referred to as the coherence. 

H. COMPUTATION OF THE PHASE ANGLE 

The average angular difference by which the cross- 

correlated points in y(t) lead their corresponding components 

of x(t) is called the cross-spectral phase angle, e (f). 

xy 

It is computed from the co- and quad-spectra by: 



Since this study deals with a time step of one day, the 
phase angles can be converted to a lag in days from: 



When either the value of the phase angle (measured in 
degrees) or the lag time (in days) is positive, then y(t) 
lags behind x(t) by that amount; if the value is negative, 
then x(t) lags behind y(t). 



e xy (f) = arctan [Q xy (f) / C xy (f)] 



(25) 



L (days ) - ^ 



(f) * 0.01745 



(26) 



f*2TT 



36 



In the computations involved in' this study there were 
8766 data points. A standard lag of 876, 10% of the total 
time series, was used. The 90% confidence interval for the 
energy densities was computed by the formula: 

r = {j- 3.7 (27) 

when N = 8716, M = 876. 

The limits were found from tables to be: 
a = 24.075 
b = 52.192 

where a and b are the values which when subtracted and 
added, respectively, to the value of the spectral energy 
density result in the lower and upper 90% confidence limits. 
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III. RESULTS 



A. REGRESSION ANALYSIS 

The regression analysis phase of this study was limited 
to the three-year time series of the NODC hydrocast data due 
to the desire to utilize the salinity data present in this 
data set alone. The period covered was from 1968 through 
1970, the same time span of data analyzed by Hansen (1973). 

At this point it is appropriate to visually inspect the 
surface temperature and salinity data sets in a determinis- 
tic approach. 

Figure 4 shows the monthly means of the sea-surface 
temperature for 1968 through 1970 compared to Robinson's 
(1971) twenty-year long-term monthly means. A characteristic 
annual cycle is evident in all three years as well as the 
long-term mean. The only major difference from the long- 
term mean is a persistent cold anomaly which is present in 
the spring of all three years. 

The year 1968 was the most anomalous, being both colder 
in March and April and warmer for the balance of the year. 
Namias (1971) attributed this to bulk heating, advection 
and mixing. The year 1969 was more normal, except for an 
anomalous cooling in the summer, attributed by Roden (1970) 
to a southward displacement of the transition zone. This 
has been supported by Wickett and Thompson (1971a, 1971b) 
whose meridional transport data. Figure 5, shows a transport 



38 




w 

<D 

3 a 

CD I 
CX. W) 

a § 

£ i 

i— t 

u r- 

ci O'' 

, -H r-H 

CO J° 



d 

o 

w 

d 

•h 

rQ 

O 



"C 
0 ) 

> 

CD 
w 

o * 

<u •§ 

5 5 

<-H O 

o r- 



CO 

d 

CD 



& 
b€ 
d 

o 

rd ^ 



CD 

U 

d 

W) 



Q 0 , ajnjoJ»diu«| 



39 



nthly Means. 




o 

m 

co 



c 

V 

V 

£ 

4 -> 

o 

pq 



u 

o 

a 

w 

C 

cd 

u 

H 



C 

o 



Jh 

V 



pH 

A 

4 -> 

§< 



nj 

4 — > 

o 

H 



in 

o 

u 

0 

ttf) 

♦ r-< 



40 



and 145 West at 30° North during 1969 and 1970. 
Positive Northward, Negative Southward (after 
Wickett and Thomson, 1971a, 1971b), 



of cooler water equatorward during April through October, 
1969. Except for the spring months, 1970 was more typical 
to the long-term mean than the preceding two years. 

Figure 6 shows the monthly means of surface salinity 
for 1968 through 1970. The most obvious fact from this 
figure is that the salinity shows evidence of having a 
semi-annual periodicity. In addition to this, there is the 
feature in 1970 of a reduction in the salinity, possibly due 
to the advection of abnormally fresh, cool subarctic water, 
indicating an equatorward shift in the subarctic transition 
zone . 

As stated previously, Hansen (1973) did much work on 
the same data series. Figure 7 is the result of his 
bivariate analysis of the monthly mean surface temperatures 
and salinities. It shows that the temperature follows four 
phases, namely the seasonal phases of summer, autumn, winter 
and spring, with autumn and spring appearing as transitional 
phases. Hansen also performed a regression test of the 
seasonal surface temperature versus the surface salinity and 
also the temperature at 50 m. The correlation coefficients 
are given in Table I. For most cases the correlations were 
poor, especially between the surface salinity and the surface 
temperature. The only good correlations occurred during the 
transition periods when both the surface temperature and 
salinity were responding to the changes of the annual cycle. 
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1968 through 1970 (from Hansen, 1973) . 



TEMPERATURE (°C) 



23 — 



22 — 



21 - 



20 - 



19 — 



18 — 



17 — 



Figure 




34.8 35.0 35.2 35.4 35.6 35.8 

SALINITY 

Monthly Temperature and Salinity Modal Points 
Resulting from the Bivariate Analysis of the 
Data, 1968 through 1970 (from Hansen, 1973). 
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Table I 



Correlation coefficients, by season, resulting from a linear 
regression analysis of surface temperature, surface salinity, 
and 50 meter temperature data, 1968 through 1970 (from 
Hansen , 1 973). 

1968 1969 1970 Mean 





Sfc 

Sal 


50 m 
T emp 


Sfc 

Sal 


50 m 
T emp 


Sfc 

Sal 


50 m 
T emp 


Sfc 

Sal 


50m 
T emp 


Winter 


0.61 


0.58 


0.61 


0.56 


0.77 


0.42 


0.66 


0.52 


Warmi ng 


0.59 


0.90 


0.77 


0.70 


-0.11 


0.67 


0.42 


0.75 


Summer 


-0.28 


0.31 


0.16 


0.40 


-0.05 


-0.29 


-0.05 


0.14 


Cooling 


0.57 


0.88 


0.84 


0.98 


-0.01 


0.96 


0.47 


0.94 


Since 


the 


surface 


salinity and 


temperature 


were 





obviously responding to different causal parameters, a 
further regression analysis was performed. The surface 
salinity was regressed against many of the classical para- 
meters associated with its fluctuations. In this instance, 
there were seven variables chosen; the vapor pressure 
difference, E-P, the heat of evaporation, the wind speed, 
evaporation, precipitation and surface pressure. No separa- 
tion into seasons was made, all computations being made as a 
three-year series of daily average observations. 

The results of the regression study are shown in Table II 
Two major observations are apparent from this study. First, 
the standard parameters often associated with a change in 
salinity, especially E-P, showed very little correlation with 
the surface salinity. However, this is quite reasonable 
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Table II 



Correlation coefficients, all parameters vs surface salinity. 



PARAMETER 


1968 


1969 


1970 


e -e 
w a 


.094 


.043 


.097 


E-P 


.045 


.151 


-.016 


QE 


.028 


.161 


-.001 


W 


-.030 


.163 


-.043 


PS 


.830 


.803 


.803 


E 


.034 


.166 


.004 


P 


-.030 


-.019 


.028 



considering the location of OWS NOVEMBER with itsyear round 
excess of evaporation over precipitation. The other major 
observation is that only one parameter is consistently 
highly correlated, i.e., the surface pressure. 

This latter fact could mean that the changes in the 
surface salinity are largely due to an advective process, 
as is supported by the fact that the other non-conservative 
terms associated with the air-sea interface have such poor 
correlation coefficients. This suggests that the changes in 
the surface salinity are related to some geostrophic process, 
most likely an advection of water from the subarctic region, 
corresponding to a southward shift of the transition zone. 
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B. SPECTRAL ANALYSIS 



The major thrust of this study was to use spectral 
analysis as a means of numerically describing the inter- 
actions at the air-sea interface at OWS NOVEMBER. Specifi- 
cally, it was to identify those meteorological parameters 
most closely associated with fluctuations in the sea-surface 
temperature. Since this latter term has a dominant annual 
cycle and most of the other parameters were compared with 
this term, the results presented here will primarily be 
those for the yearly cycle. The logic of this is further 
demonstrated by comparing the percentages of the variance 
explained by the yearly cycle for the different terms. The 
percentages shown in Table III were calculated from the 
energy-density functions which were produced for each time 
series. This function is a measure of the variance contained 
at that frequency. To compute the variance contributed by 
the yearly cycle a bandwidth was specified by those frequencies 
for which the value of the energy-density function decreased 
to less than half the value at the yearly cycle. 

From Table III it can be seen that 86% of the total 
variance of the sea-surface temperature is explained by the 
yearly cycle alone. For this reason, even though the 
spectral results were computed for frequencies as short as 
the monthly cycle, the emphasis of this study was placed on 
the yearly cycle. Also, since the major areas of interest 
were the correlations with the heat budget terms, each of 
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these is discussed describing its influence on the sea- 
surface temperature. In addition, the contributing effects 
of the parameters involved in the computation of that heat 
flux term will be discussed. 



Table III 

Contribution of the yearly cycle to the total variance explained 



Parameter 

SST 

T a 

T w - T a 
e w 



e -e 
w a 

E (evaporation) 

W 

w 2 

Cl (cloud cover) 

PS 

QS 

QB 

QE 

QC 

QN 



% of Total Vari ance 
86 
80 
20 
86 
46 
14 
34 
8 
9 

7 

8 
70 
87 
11 
17 
42 
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The first heat budget term to be examined is the solar 
radiation, QS , which like the surface temperature exhibits 
a strong annual periodicity. From Tabe III, it can be seen 
that 70 % of the total variance of this term is explained by 
the annual cycle. The term which primarily controls the 
variation of QS is the cloud cover. Cl, and in Table IV the 
spectral results for the yearly cycle are given for the 
interactions of the SST and Cl with QS. Along with a 
correlation of the cloud cover with the SST, the complete 
spectral results are shown in Figures 8, 9, and 10. Figures 
A-l , A-2, and A-3 of Appendix A show the respective auto- 
correlation functions and Figures B-l and B-2 of Appendix B 
illustrate 15-year segments of the 24-year time series of 
the SST and QS. 

As would be expected. Table IV shows that the solar 
radiation is strongly coherent with the SST at the yearly 
cycle. However, there is a lag of almost 60 days, which 
cannot be explained by the cloud cover correlation with its 
low coherence and therefore dubious lag time. Also, there 
is a strong coherence at the semi-annual cycle, which appears 
to reflect the effect of the cloud cover as seen in Figure 10. 



Table IV 



Spectral results for the solar radiation, 
SST and cloud cover at the yearly cycle. 



QS vs SST 



Parameter 

Codes 



Coherence 



.94567 



Phase 
(degrees ) 



- 60.0 



Lag 

(days ) 
-58.5 



QS vs Cl 



. 85236 



126.2 



122.0 
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OHCRENCE LOG SPECTRAL ENERGY DENSITY 



o 




TIME(DAYS) 
f (DAYS' 1 ) 



TIME (DAYS) 
f (DAYS' 1 ) 




Figure 8. Spectral density, coherence and phase for 
the SST vs the solar radiation. 
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PH9SE COHERENCE L3G SPECTR9L ENEPGV DENSIT V 

-ICQ Q ]CQ O U 1 O -Z -0 z 



F -I 

(Days ) 


TIME 
(Days ) 


0.007 = 


143 


0.013 = 


77 


0.020 = 


37 


0.034 = 


29 



SST 




.000 



0.007 



0.013 



0.020 



0.027 



0.034 



TIME(DAYS) 
f (DAYS*’ ) 





350 175 143 77 50 37 29 

I 1 H 1 1 1 1 

000 0.007 0.013 0.020 0.027 0.034 



TIME(DAYS) 
f (DAYS' 1 ) 




Figure 9. Spectral density, coherence and phase for 
the SST vs cloud cover. 
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77 

| 


50 

1 


37 

1 


29 

- —1 


TIME (DAYS) 


1 

0.007 


1 

0.013 


0.020 


1 

0.027 


1 

0.034 


f (DAYS' 1 ) 




TIME (DAYS) 
f (DAYS*’) 




Figure 10. Spectral density, coherence and phase for 
the solar radiation vs cloud cover. 
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This rather large lag is equally evident in Table V 
which shows the months in which SST and the heat transfer 
terms achieve their maximum and minimum values. Along with 
Figures B-l and B-2, it can be seen that the surface 
temperature is a maximum during September and at a minimum 
during March. However, the solar radiation peaks in August 
and is a minimum in January. This is approximately the 58 
day lag shown in Table IV. This lag can be explained by the 
fact that the solar radiation is absorbed and transferred 
throughout the upper mixed layers of the water column. Thus, 
a lag of two months exists before the total effects of the 
solar radiation are represented by a warming or cooling of 
the surface temperature. 



Table V 

Months of maximum and minimum value for the 
surface temperature and the heat budget term. 



Parameter 



Mi ni mum 



Maxi mum 



SST 


March 


September 


QS 


January 


August/September 


QE 


May 


Novembe r 


QB 


July/January 


March/November 


QC 


August 


March 


QN 


December 


June 
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In Figures 8, 9, and 10 there are also noticeable peaks 
at the six month cycle and three month cycle. These results 
are thought to be, for the most part, artifacts of the yearly 
cycle and are relatively small when compared to the yearly 
value; yet, as mentioned before, several show evidence of 
being real . 

The major result of this portion of the study is that 
the use of the new formula for QS results in significantly, 
higher estimates of this term than previously calculated. 
Fluctuations in solar radiation precede those of the surface 
temperature by approximately two months and thus the solar 
radiation is a good predictor of the surface temperature, 
as would be expected. 

The back radiation, QB, represents the loss of heat 
energy as long-wave radiation. It is strongly dependent 
upon the surface temperature and the vapor pressure of the 
atmosphere. Table VI summarizes the results for the yearly 
cycle of QB and associated terms with the surface temperature. 

Table VI 

Spectral results for the back radiation and 
related parameters at the yearly cycle. 



QB vs SST 


.59755 


121.8 


119.0 


T vs SST 

a 


.98800 


- 6.0 


- 5.9 


T.-T vs SST 
w a 


.67194 


-49.1 


-47.9 


e, vs SST 

a 


.98815 


- 6.4 


- 6.2 


QB vs Cl 


.37565 


-82.9 


-78.1 
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From Table III it was found that about 87% of the total 
variance of the back radiation was contained in the yearly 
cycle; yet, Table VI shows only a slight coherence between 
the SST and QB for the yearly cycle and a large dubious lag 
value. This can be explained by examining the interactions 
involved. The back radiation term decreases as the SST 
increases due to the effect of the humidity in the atmosphere 
directly above the surface. Thus, over the yearly cycle, 

QB would ideally be 180 degrees out of phase with the surface 
temperature. However, the back radiation is also decreased 
by an increase in the cloud coverage. Thus, the condition 
illustrated by Figure B-3 and Table V may be seen where QB 
is shown to have two maxima and minima on roughly a semi- 
annual basis. From Figure A-4 it may be seen that this 
feature is superimposed over a larger annual cycle which 
explains the high percentage of the variance contained by the 
yearly cycle. 

The complete spectral results for the above correlations 
are given in Figures 11 through 15. From Figure 11 it is 
seen that QB has two peaks, one at the yearly cycle and 
another lesser one at one-half the yearly cycle. Of the 
related parameters, only the air temperature and the vapor 
pressure of the air are highly correlated and then at only 
the yearly cycle. The air-sea temperature difference is only 
mildly correlated at the yearly cycle and its associated lag 
is highly suspect for that reason. The air-sea temperature 
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PHASE COHERENCE LOG SPECTRAL ENERC V DENSITY 




O V. 



350 175 143 77 50 37 29 

I 1 H I 1 1 1 

0.000 0.007 0.013 0.020 0.027 0.034 




Figure 11. Spectral density, coherence and phase for 
the SST vs the back radiation. 
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PHRSE COHERENCE L0G SPECTRRL ENERGV DENSITY 





Figure 12. Spectral density, coherence and phase for 
the SST vs air temperature. 
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Figure 13. Spectral density, coherence and phase for 

the SST vs the air-sea temperature difference. 
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Figure 14. Spectral density, coherence and phase for the 
SST vs the vapor pressure of the air. 
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Figure 15. Spectral density, coherence and phase for the 
back radiation vs cloud cover. 



59 



difference demonstrates a coherence at the 70-day cycle. 

Also the back radiation and cloud cover show evidence of 
this possible seasonal coherence. 

The major conclusion to be drawn from this part of the 
study is that the back radiation is semi-annual with an 
underlying yearly cycle. As a predictor for the surface 
temperature, the back radiation is poorly associated with 
fluctuations in SST. 

The next heat budget term examined was the heat of 
evaporation, QE. Only 11% of the total variance of QE is 
explained by the yearly cycle which in itself means that QE 
is a poor predictor of the surface temperature. This is 
verified by Table VII which shows only a moderate coherence 
for QE and the SST at the yearly cycle. This lessens the 
credibility of the large lag shown; also from Table V it can 
be seen that the actual lag is more of the nature of two 
months. Additionally the lag is positive which means that 
fluctuations of QE follow those of the SST. 



Table VII 

Spectral results for the heat of evaporation 
and related parameters at the yearly cycle. 



Parameter 

Codes 


Coherence 


Phase 
( degrees ) 


Lag 

(days ) 


QE vs SST 


.84783 


97.5 


95.0 


e -e vs SST 
w a 


.86425 


22.8 


22.3 


W vs SST 


.95593 


138.1 


135.4 


W 2 vs SST 


.95512 


135.4 


132.0 


E vs SST 


.81722 


76.9 


74.9 
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Whereas QE is only weakly correlated with the SST at 
the annual cycle, the wind .speed and wind speed squared are 
highly correlated to the SST with large positive lags 
(approximately 4 months). These would seem to explain much 
of the lag in QE when they are coupled with the vapor 
pressure difference. QE is a transfer mechanism of the form: 

QE = - A e || (28) 

where : 

A g is the eddy coefficient for water vapor 

9 f 

— is the humidity gradient above the sea surface. 

Both of these are highly dependent upon the wind speed. 

Even though the vapor pressure difference and evaporation 
are weakly correlated to the SST, they combine with the wind 
speed to produce the low coherence and large lag for QE . 

Figures 16 through 20 give the complete spectral results 
for the correlations shown. Being computed from the surface 
temperature, e w and SST are naturally perfectly correlated 
and is not included. The yearly cycle dominates in all 
except for QE, the vapor pressure difference and the evapora- 
tion, as would be expected. Both the wind speed and wind 
speed squared show a definite coherence at about the 80-day 
cycle. The autocorrelation functions are found in Figures 
A-9 through A - 1 3 and a 15-year plot of QE and its anomaly 
is shown in Fi gure B-4 . 
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Figure 16 . Spectral density, coherence and phase for the 
SST vs the heat of evaporation. 
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Figure 18. Spectral density, coherence and phase for the 
SST vs the wind speed. 
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Figure 19. Spectral density, coherence and phase for the 
SST vs the wind speed squared. 
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Figure 20. Spectral density, coherence and phase for the 
SST vs the evaporation. 
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The heat of conduction, QC, is the smallest of the heat 
flux terms and only 17% of its total variance is explained 
by the yearly cycle. From Table VIII it can be seen that 
this term also has an unusually large lag, approximately 
130 days, and the lowest coherence of any of the heat budget 
terms with the surface temperature. From Figure 21 it can 
be seen that QC is nearly of the form of white noise with 
no distinct periodicities. This is confirmed by Figure B-5 
which shows that QC is constantly small and contains few 
anomalies of note. 

Like QE , QC is a flux term of the form 

QC = -A h fl (29) 

wiie re 

A^ is the eddy coefficient of heat conduction, and 

aT 

t— is the temperature gradient of the air overlying 
the sea surface. 

For an area of the nature of OWS NOVEMBER, dominated by 
fairly constant temperature patterns for both the air and 
the ocean, it can be seen that this term would be expected 
to be fairly constant. It does contain quite a deal of 
energy in cycles other than the annual and shows evidence 
of good coherence at the semi-annual, seasonal and monthly 
peri odi ci ti es . 
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TABLE VIII. 



Spectral results for the convective 
heat vs the SST at the yearly cycle 



Pa rameter 
Codes 



Coherence 



Phase 
(degrees ) 



Lag 
(days ) 

130.7 



QC vs SST 



. 75694 



134.1 



The final heat budget term, the net heat transfer, QN , 
is perhaps the most interesting of the heat equation terms, 
representing the gain or loss of heat by the water column. 

To this point, the advective and local heat exchange terms 
have been ignored. However, since the values for the net 
radiation computed from Equation (1) did not result in zero 
values, then it is evident that these terms have to be 
considered. If QN is then viewed as a combination of the 
advective heat term and the local generation of heat it takes 
on the form: 



where : 

Q v is the advection of heat by currents, and 

Q q is the local heat generated. 

0 

Over the annual cycle Q_ = 0; therefore, QN must be 
representative of only the advection term. This is supported 
by the fact that for the entire 24-year period QN was found 
to have a slope of only 0.003 which signifies that no long 
term generation of heat has occurred. 



ON = Q v - Q 0 
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From Tables III and IX it can be seen that the net heat 
exchange is dominated by the annual cycle and is highly 
coherent with the surface temperature at this period. 
Forty-two percent of the total variance is explained by the 
yearly cycle, which allows for some seasonal effects and 
Figure 22 shows some evidence of high coherence at the 
semi-annual and seasonal cycles. The lag is confirmed by 
both Table V and Figure B-6. Much of this lag is undoubtedly 
due to the effect of the 58-day lag of the solar radiation 
and the balance by the other heat budget terms. In addition, 
the lag of QN and the surface pressure, considered here 
because of the previously noted high correlation of pressure 
with the surface salinity, may also contribute to the lag of 
QN and SST. 



TABLE IX 

Spectral results for the net radiation 
and related parameters at the yearly cycle 



Parameter 

Codes 


Coherence 

Squared 


Phase 
(degrees ) 


Lag 

(days ) 


QN vs SST 


.97954 


- 76.7 


- 74.8 


PS vs SST 


. 79993 


-120.6 


-118.6 


QN vs PS 


.79250 


- 46.9 


- 45.3 
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Figure 22. Spectral density, coherence and phase for the 
SST vs the net radiation. 
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Figures 22 through 24 give the complete spectral 
results for the net heat exchange, the surface temperature, 
and the surface pressure correlated versus each other. In 
all three instances the yearly peak dominates. 

The major result of this segment is that the new solar 
radiation formula results in larger positive values of QN . 
Secondly, for the yearly cycle QN is a measure of heat 
advection and at OWS NOVEMBER this advection is of nearly 
equal importance as the solar radiation as an influence on 
the surface temperature. The negative lag makes QN a fairly 
good predictor of the surface temperature, even though much 
of the variance of the net heat exchange is contained in the 
seasonal cycles. The surface pressure seems to be less 
correlated with the surface temperature than it is with the 
salinity. 

Very little has been mentioned about the autocorrela- 
tions; however, White (1974) has used this function in the 
study of SST and baroclinic transport anomalies.. He found 
an apparent 5-year periodicity in the SST anomalies in the 
vicinity of OWS NOVEMBER. This cycle was attributed to 
deeper anomalous features in the baroclinic structure. 
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Figure 23. Spectral density, coherence and phase for the 
SST vs the surface pressure. 
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Figure 24. Spectral density, coherence and phase for the 
net radiation vs the surface pressure. 
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IV. CONCLUSIONS 



The aim of this study was to use the tools of spectral 
analysis and other statistical approaches to describe the 
interactions of the air-sea interface at OWS NOVEMBER, 
specifically as to finding those parameters most closely 
associated with fluctuations of the daily average surface 
temperature and surface salinity. 

The results of this research showed that: 

1. The surface salinity from 1968 through 1970 was 
shown through regression analysis to be poorly correlated 
with the parameters traditionally associated with changes 
in the observed salinity. The evaporation, evaporation 
minus precipitation, and precipitation were found to have 
little correlation with fluctuations in surface salinity. 
These terms are non-conservative terms associated with the 
air-sea interface and their low correlation coefficients 
would appear to negate their importance in controlling the 
surface salinity at this location in the Pacific. 

2. The surface pressure was found to be the most 
strongly correlated of all the parameters regressed against 
the surface salinity. It was therefore hypothesized that 
this represents the presence of a movement of a baroclinic 
system of cooler water from the subarctic region into the 
vicinity of OWS NOVEMBER, corresponding to a shift southward 
of the transition zone. 
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3. For the yearly cycle, it was shown that the net 
radiation is a measure of the advection present and that 
this advection is of nearly equal magnitude as the solar 
radiation. 

4. The new formula for the solar radiation leads to 
values for this term approximately 20 to 30 percent higher 
than previously computed. In turn, this leads to much 
larger positive values of the net heat exchange term and 
larger positive anomalies for both terms. 

5. From the spectral analyses it was found that all 
the 24-year time series were dominated by the yearly cycle; 
except for the back radiation, cloud cover, wind speed and 
wind speed squared series which showed evidence of semi- 
annual and seasonal cycles as well. 

6. As predictors of the surface temperature, the 
solar radiation and advection in the form of the net heat 
exchange term were the only terms of suitable phase and 
coherence at the yearly cycle to be valid. The o.ther heat 
transfer terms seems to have less effect upon fluctuations 
of the surface temperature. Several of the observed para- 
meters, such as the air temperature, were also highly 
coherent and, on the whole, had smaller lags associated with 
them . 
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APPENDIX A 



THE AUTOCORRELATION FUNCTIONS 



Figure Page 

A-l Autocorrelation curves for the SST 

and the solar radiation, QS 78 

A-2 Autocorrelation curves for the SST 

and cloud coverage 79 

A-3 Autocorrel ati on curves for the cloud 

coverage and the solar radiation 80 

A-4 Au tocorrel ati on curves for the SST 

and the back radiation 81 

A-5 Autocorrelation curves for the SST 

and the air temperature 82 

A-6 Autocorrel ati on curves for the SST 

and the air-sea temperature difference 83 

A-7 Autocorrel ati on curves for the SST 

and the vapor pressure of the air 84 

A-8 Autocorrelation curves for the 

back radiation and cloud cover 85 

A-9 Autocorrel ati on curves for the SST 

and the heat of evaporation 86 

A-10 Autocorrel ati on curves for the SST 

and the vapor pressure difference 87 

A - 1 1 Autocorrel ati on curves for the SST 

and the wind speed 88 

A - 1 2 Autocorrel ati on curves for the SST 

and the wind speed squared 89 

A - 1 3 Autocorrel ati on curves for the SST 

and evaporation 90 

A - 1 4 Autocorrel ati on curves for the SST 

and the convective heat 91 

A - 1 5 Autocorrelation curves for the SST 

and the net radiation 92 

A - 1 6 Autocorre 1 ati on curves for the SST 

and surface pressure 93 

A - 1 7 Autocorrel ati on curves for the net 

radiation and the surface pressure 94 
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Figure A-l. Autocorrel ati on curves for the SST and the solar 
radiation, QS (abscissa is given in hundreds 
of days ) . 
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Figure A-2. Autocorrelation curves for the SST and cloud 

coverage (abscissa is given in hundreds of days). 
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Figure A-3. Autocorrelation curves for the cloud coverage 
and the solar radiation (abscissa is given in 
hundreds of days ) . 
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Figure A-4. Autocorrelation curves for the SST and the back 

radiation (abscissa is given in hundreds of days). 
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Figure A-5. Autocorrelation curves for the SST and the air 

temperature (abscissa is given in hundreds of days). 
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Figure A-6. Autocorrelation curves for the SST and the 
air-sea temperature difference (abscissa is 
given in hundreds of days.). 
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Figure A-7. Autocorrelation curves for the SST and the vapor 
pressure of the air (abscissa is given in 
hundreds of days ) . 
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Fiqure A- 8. Autocorrelation curves for the back, radiation 
and cloud cover (abscissa is given in hundreds 
of days). 
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Figure A-9. Autocorrelation curves for the SST and the heat 
of evaporation (abscissa is given in hundreds 
of days). 







Figure A - 1 0 . Autocorrelation curves for the SST and the vapor 
pressure difference (abscissa is given in 
hundreds of days ) . 
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Figure A - 1 1 . Autocorrelation curves for the SST and the 
wind speed (abscissa is given in hundreds 
of days ) . 
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Figure A - 1 2 . Autocorrelation curves for the SST and the 
wind speed squared (abscissa is given in 
hundreds of days). 
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Figure A-13. Autocorrelation curves for the SST and 

evaporation (abscissa is given in hundreds 
of days). 
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Figure A- 1 4 . Autocorrel ati on curves for the SST and the 

convective heat (abscissa is given in hundreds 
of days ) . 
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Figure A - 1 5 . Autocorrelation curves for the SST and the net 
radiation (abscissa is given in hundreds of 
days ) . 
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Figure A-16. Autocorrelation curves for the SST and surface 

pressure (abscissa is given in hundreds of days). 
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Figure A-17. Autocorrelation curves for the net radiation 
and the surface pressure (abscissa is given 
in hundreds of days ) . 
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APPENDIX B 



FIFTEEN-YEAR TIME SERIES SEGMENTS 

Figure Page 

B-l Fifteen-year time series of the 

surface temperature and its anomaly 96 

B-2 Fifteen-year time series of the 

solar radiation and its anomaly 9 7 

B-3 Fifteen-year time series of the 

back radiation and its anomaly 98 

B-4 Fifteen-year time series of the 

heat of evaporation and its anomaly 99 

B-5 Fifteen-year time series of the 

convective heat and its anomaly 100 

B-6 Fifteen-year time series of the 

net radiation and its anomaly 101 
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Figure B-2. Fifteen-year time series of the solar 
radiation and its anomaly (cal cm -2 d 
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Figure B-4. Fifteen-year time series of the heat 

of evaporation and its anomaly (cal cm 
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